Two-proton removal from 44 S and the structure of 42 Si 
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Newly-published 42 Si gamma-ray spectra and a final-state-inclusive 42 Si production cross section 
value, obtained in a higher-statistics intermediate-energy two-proton removal experiment from 44 S, 
are considered in terms of the final-state-exclusive cross sections computed using proposed shell- 
model effective interactions for nuclei near N = 28. Specifically, we give cross section predictions 
when using the two nucleon amplitudes of the two-proton overlaps ( 42 Si( J") | 44 S) computed using 
the newly-proposed SDPF-mu shell-model Hamiltonian. We show that these partial cross sections 
or their longitudinal momentum distributions should enable a less-tentative interpretation of the 
measured gamma-ray spectra and provide a more quantitative assessment of proposed shell-model 
Hamiltonians in this interesting and challenging region of the chart of nuclides. 

PACS numbers: 24.50.+g, 25.60.-t, 25.70.-z, 27.40.+Z 



I. INTRODUCTION 

Measurements and theoretical predictions of the prop- 
erties of the neutron-rich N = 28 isotones are subjects 
of considerable interest and recent studies, in an attempt 
to understand quantitatively the role of cross-shell inter- 
actions on the erosion of the TV = 28 sub-shell gap - 
particularly for nuclei in the vicinity of 42 Si. This in- 
terest, the major issues involved, and much of the key 
literature are summarised concisely in the recent Letter 
by Takeuchi et al. pj, that also reports a state-of-the- 
art, higher-statistics, in-flight 7-ray spectroscopy mea- 
surement of final states of 42 Si. The experiment, using 
a fast 44 S(— 2p) two-proton removal reaction, was per- 
formed at the RIBF, RIKEN, at 210 MeV per nucleon. 
Accessing this region near Z = 14 is a considerable ex- 
perimental challenge. The region also provides new and 
severe challenges for theoretical calculations. 

The spectrum for 42 Si deduced in Ref. [l| was com- 
pared with the level-scheme from shell-model calcula- 
tions computed from a newly-constructed sd — pf cross- 
shell Hamiltonian, named SDPF-MU @]. This very re- 
cent theoretical study includes and shows very clearly 
the importance of the tensor-force component in the 
cross-shell sd — pf proton-neutron interaction for repro- 
ducing the measured 2+- and 4+-state energies and the 
B(E2,0 + — > 2+) values in the silicon and sulfur isotopes 
as N approaches 28 - including the energy of the (tenta- 
tively) assigned 42 Si 4+ state in reference 

The analysis of Ref. [l| tentatively assigned the ob- 
served 1431(11) keV 7-rays, also seen in a coincidence 
gate with the 742(8) keV 7-ray, previously identified as 
due to the 2^ — > + (g.s.) transition 0], to the decay of 
the 42 Si(4^) state located at 2173(14) keV. The Af state 
is at 2270 keV according to the SDPF-MU calculation. 
This and the other SDPF-mu bound 42 Si( J T ) final states 
below the (somewhat uncertain) first neutron threshold, 
of 3.6(6) MeV 0], are collected in the first two columns 



of Table [U The additional states up to 4.2 MeV are also 
shown. These shell-model calculations include the full 
sd-shell (for protons) and the full p/-shell (for neutrons) 
and used the code NuShellX There are candidate 
sdpf-mu shell-model states for each of the decay 7-rays 
observed in the inclusive spectrum of [l[ since the 4+, 
O2 and 0^ shell-model states all decay through the first 
excited 2J~ state. 

In contrast, a shell- model calculation that uses instead 
the SDPF-u-Si variant of the SDPF-u interaction, designed 
for Oftw-truncated calculations of the neutron-rich silicon 
isotopes in the same model space as is used here @ , pro- 
duces a more fragmented spectrum with 23 sub-threshold 
bound J v = + , . . . , 6 + , 42 Si final states, 21 of which can 
be populated directly by the removal of two sd-shcll pro- 
tons. These states, up to a threshold energy of 3.6 MeV, 
arc shown in Table \H\ Two 1 + states, at 2985 and 3380 
keV, that are found to have negligible direct 2p-removal 
cross sections are omitted from the Table, as are the pre- 
dicted bound 5+ (3288 keV) and 6+ (3472 keV) states 
that cannot be populated directly. There are a further 9 
excited states in the interval from 3.6 to 4.2 MeV. 

In the following we calculate the cross sections for pop- 
ulation of the individual 42 Si( J T ) final states in each case, 
resulting from the two-proton removal reaction on a car- 
bon target at 210 MeV per nucleon, the reaction of [l|. 



II. CROSS SECTION CALCULATIONS 

The removal of two like nucleons of the deficient species 
from highly N : Z asymmetric nuclei has: (a) been shown 
to proceed directly [7j and to provide spectroscopic in- 
formation [H, and (b) has been assessed quantitatively 
using several 2p- and 2n-removal reaction test-cases on 
sd-shell nuclei @, Q . We follow the approach described 
in detail in these latter references. The structure de- 
tails from the shell-model calculations enter the removal 
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reaction calculations through the two-proton overlaps, 



i.e. 
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1,2) = (1,2, 42 Si( JM) | 44 S(0+)). Specifically, 
for each final-state, E x , J, they generate the set of two 
nucleon amplitudes (TNAs) Cj^- 2 for each active two- 
proton configuration in the assumed shell-model space 
where, for a spin + projectile, 



J+M 



I J [<Ah (1) ® 4>h (2)] J-M • 



Given these TNAs, the reaction is described as a sud- 
den, direct removal process in which the removed pro- 
tons interact both elastically and inelastically and reac- 
tion residues only elastically with the light nuclear tar- 
get, here carbon. These absorptive two-body interactions 
arc described by their elastic S-matrices, computed using 
the eikonal model from the complex residue- and proton- 
target optical potentials. These potentials, at 210 MeV 
per nucleon, were computed using the double- and single- 
folding models, respectively, assuming a zero-range effec- 
tive nucleon- nucleon interaction based on the free pn and 
pp cross sections. The point nucleon density of the car- 
bon nuclei was assumed to be of Gaussian form with a 
root mean squared (rms) radius of 2.32 fm. That of the 
42 Si residues is discussed below. 

As has been detailed for a general case in Ref. [l(| , the 
following nuclear model parameters were used, (a) The 
point neutron and proton densities of the 42 Si residues 
were obtained from spherical Hartrec-Fock calculations 
using a Skyrme force (the SkX interaction [H[). The 
removed-proton se?-shell radial form factors were de- 
scribed as normalized eigenstates of a Woods-Saxon plus 
spin-orbit potential well whose depths were constrained 
by the evaluated 2p separation energy S2 P = 40.2(5) 



TABLE I: Spins, parities and excitation energies of all shell- 
model states below the neutron separation energy of 42 Si, 
3.6(6) MeV, computed using the SDPF-mu Hamiltonian of Ref. 

The additional states up to 4.2 MeV excitation energy are 
also shown. The two-proton removal partial cross sections to 
each final state and the inclusive cross section for states up 
to a 3.6 MeV threshold are also shown, calculated using the 
SDPF-MU TNAs . 



-Esdpf-mu (keV) 



oth(-2p) (mb) 



or 
2?- 

0+ 

oi 

2 2 + 

Inclusive cr t h( — 2p): 

2 3 + 
3T 
4T 
2* 






0.186 


821 


0.031 


2271 


0.030 


2573 


0.103 


3273 


0.011 


3525 


0.012 




0.37 


3844 


0.005 


3899 


0.007 


4080 


0.010 


4090 


0.007 



TABLE II: Spins, parities and excitation energies of all 42 Si 
shell-model states below a 3.6 MeV first neutron threshold 
energy having non- negligible direct two-proton removal reac- 
tion yields. The 42 Si states were computed using the SDPF- 
U-Sl interaction For the two-proton removal TNA and 
partial cross section calculations, the wave function for 44 S 
was computed using the SDPF-u interaction. 
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S P = (S2p+E x )/2. We use excitation energies E x 



i.e. 

E SM 

as given by the shell-model. Fixed diffuseness (0.7 fm) 
and spin-orbit interaction strength (6.0 MeV) parame- 
ters were assumed throughout. The radius parameters 
of these Woods-Saxon binding potentials, r , were de- 
termined from a fit to the rms radius of each proton 
single-particle orbital of interest, as given by the HF- 
SkX calculation. These rms radii play a significant role 
in the determination of the removal cross sections, since 
they dictate the spatial extent of the two-proton position 
probabilities relative to the 42 Si residue densities, and 
these two sizes must be specified consistently. Hence our 
usual practice of using the same theoretical (HF) model 
to constrain the proton radii and the densities used in 
the evaluation of the interaction S'-matrices. 

The 42 Si residue longitudinal momentum distributions 
offer an additional signal of the spins of the populated 
and 7-decaying final-states, as has been discussed for- 
mally [HI, EH and also applied for spectroscopy; for ex- 
ample, in identifying the excited 2t and 4^ states in 44 S 
following the 46 Ar(-2p) reaction [lj|- In the following 
section we will discuss the exclusive momentum distribu- 
tions of selected SDPF-mu shell-model final states. Those 



3 



calculations use the same input parameters as have been 
detailed above. 



III. RESULTS 

The individual final-state partial cross sections and 
longitudinal momentum distributions were calculated fol- 
lowing the formalism of Refs. [1,0] and [HI, EH, respec- 
tively, from the np-form shell-model TNAs computed us- 
ing NuShellX [|. 



A. Exclusive cross sections 

The calculated final-state-cxclusive two-proton re- 
moval cross sections for SDPF-mu are shown in Table U 
that shows the strongest populations are of the ground- 
state and the first excited 0^ state. The calculated 
inclusive removal cross section is 0.37 mb, to be com- 
pared with the reported value of 0.15(2) mb So, 
the ratio of the experimental to the theoretical yields, 
R B (2N) = cr cxp /(7 th (iiicl.) = 0.41(6), is consistent with 
the systematics, R S (2N) « 0.5, previously observed for 
a number of sd-shell test-case nuclei @ . 

The corresponding situation, using the TNAs from the 
wave functions for 44 S computed using SDPF-u and for 
42 Si computed using SDPF-u-Si, are shown in Table fll] 
The integrated cross section to the (21) SDPF-u-Si fi- 
nal states is 0.69 mb. The cross section to the lowest 
6 states, the 0^~ 2 3: 2^~ 2 and 4^~, in common with the 
SDPF-mu bound states (though reordered), is 0.41 mb. 
These states all lie below 1.8 MeV when using sdpf-u-si, 
indicating that significant additional strength is present 
at lower excitation energies and that is distributed over a 
large number of the additional bound final states. A very 
similar integrated cross section (of 0.62 mb) and pattern 
of final-state populations arises if instead one computes 
both the wave functions for 44 S and 42 Si using SDPF- 
U-Si. These calculated cross-section distributions show 
little correspondence with the latest experimental obser- 
vations [l|. 

The predicted relative final-state yields of Table Q] on 
the other hand suggest clear fingerprints for a more de- 
tailed assessment of the SDPF-mu wave functions from 
the measured 7-ray yields, predicting (i) that 50% of the 
inclusive cross section will pass through the 2^~ state, 
and (b) specific ratios for the proposed (4+)-state decay 
strength relative to those for the 2+-state and the pre- 
dicted excited O^-state decays. 

We note that prediction (i) above is consistent with 
the observation from the earlier lower energy and lower 
statistics 44 S(— 2p) measurement made at GANIL [H, 
which reported that 44(10)% of the reaction events were 
accompanied by a 2f to ground state transition 7-ray - 
there measured to be 770(19) keV in reasonable agree- 
ment with the new data of ll|. 
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FIG. 1: (Color online) Calculated exclusive longitudinal mo- 
mentum distributions in the projectile rest frame of 42 Si 
residues produced in the (solid curve) and 4^ (dashed 
curve) shell-model final states of the SDPF-mu Hamiltonian 
0- These calculated momentum distributions are normal- 
ized to unity at their peaks. Their integrated cross sections 
were shown in Tabled] 

B. Exclusive momentum distributions 

Further to the fingerprints offered by these relative 
partial cross sections, in Figure Q] we now show the 
computed exclusive longitudinal momentum distribu- 
tions d<Jth(J 7T )/dk c , in the projectile rest frame, for 42 Si 
residues produced in the 4^~ (dashed curve) and 0^ (solid 
curve) SDPF-mu shell-model final states. The 4f state, 
that is dominated by the ((I3/2 ® ^5/2)4 two-proton con- 
figuration with a dominant L — 4 total orbital angular 
momentum component, is particularly broad. This con- 
trasts markedly with the 0j distribution. These large 
differential widths suggest that one could confirm the ex- 
istence of these states from their momentum distributions 
in measurements with quite limited statistics. 



C. Inclusive momentum distributions 

Although less specific, the inclusive momentum distri- 
butions for (a) reactions to all final states, and (b) for all 
events with a coincident 2^-state 7-ray, both involving 
increased statistics, would also probe the SDPF-MU pre- 
dictions. The predicted distributions, obtained by sum- 
ming the appropriate exclusive distributions, are shown 
by the solid and dashed lines, respectively, in Figure [5] 
The predicted inclusive cross section, shown in Table U to 
be dominated by + transitions (81%), is predicted to be 
narrow, analogous to the solid, 0^-state curve shown in 
Fig. [TJ The 2+-state inclusive events are predicted to be 
a 60 : 40% admixture of 0+ : (2+ + 4+) transitions. This 
distribution is seen to be broadened somewhat by the 
now 40% component from the 2 + - and 4 + -state distribu- 
tions. These two distributions thus provide a valuable, 
but less specific profile of the predicted SDPF-mu bound 
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FIG. 2: (Color online) Calculated inclusive longitudinal mo- 
mentum distributions in the projectile rest frame for (i) reac- 
tions to all final states (solid curve) and (ii) for all events with 
a coincident 2^-state 7-ray (dashed curve). Calculations use 
the SDPF-mu shell-model TNAs. These momentum distribu- 
tions are normalized to unity at their peaks. 

states yields, in particular of the + -states dominance 
and, in the case of the coincident 2+-state events, the 
relative importance of excited + -states of the SDPF-MU 
Hamiltonian. 



IV. SUMMARY COMMENTS 

We have presented an analysis of the partial cross sec- 
tions for the 44 S(— 2p) two-proton removal reaction at 
210 MeV per nucleoli, populating the bound final states 
of 42 Si predicted by the newly-developed SDPF-mu cross- 
shell Hamiltonian of Ref. ■ The calculations make spe- 
cific predictions for the relative populations of the final 



states, as were shown in Tabic HI that should be tested 
against deduced 7-ray yields - as should be available from 
the data of 

This could allow a considerable strengthening of the 
tentative assignment of a 2173(4) keV 4+ excited state 
in 42 Si made in [l[ and a more detailed assessment of 
the predictions of the SDPF-mu Hamiltonian than can 
provided by, for example, the 2^~ and 4^ level energies. 
This includes the prediction of a strong population of 
an excited state near 2.5 MeV. The SDPF-MU-based 
final-state-exclusive reaction yields presented are consis- 
tent with the earlier experimental observation 0] that 
44(10)% of 44 S(— 2p) reaction events were accompanied 
by a 42 Si 2+ to ground-state decay 7-ray. 

We also show that the computed longitudinal momen- 
tum distributions predict clear additional fingerprints of 
the populated final-states, with large differential widths 
calculated for 42 Si residues produced in the 4+ and Oj" 
states, states that are predicted by the SDPF-mu Hamil- 
tonian to lie at quite similar excitation energies. The 
inclusive momentum distributions are shown to be nar- 
row reflecting the dominance of + final-state transitions 
predicted by the SDPF-MU Hamiltonian. 

In closing we note that complementary 42 Si final-state 
data from the 43 P(— lp) one-proton removal reaction 
would also contribute to this 42 Si spectroscopy discus- 
sion, since sd-shell proton removal from the 43 P(i + ) 
ground-state cannot directly populate proposed J 11 = 4+ , 
42 Si final states. 
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